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1.  Abstract 


Progress  has  been  made  on  understanding  the  complex  behavior  of  physical  processes  described 
by  nonlinear  ordinary  and  partial  differential  equations  through  the  use  of  singular  perturbation 
methods.  Modulation  equations  for  the  amplitude  and  phase  of  dissipatively  perturbed  strongly 
nonlinear  oscillators  and  traveling  waves  have  been  derived  from  the  action  equation  using  the  usual 
method  of  multiple  scales.  Equivalent  results  have  been  obtained  using  the  method  of  averaging 
developed  for  the  first  time  for  a  nonlinear  partial  differential  equation,  the  Klein-Gordon  equation, 
describing  dispersive  waves.  In  another  study,  Whitham’s  averaged  Lagrangian  principle  has  been 
generalized  to  account  for  arbitrary  perturbations  of  the  initial  conditions.  In  other  work,  Bourland 
and  Haberman  analyzed  the  slow  crossing  of  an  unperturbed  homoclinic  orbit  (separatrix)  for 
dynamical  systems.  Solutions  in  the  neighborhood  of  the  separatrix  are  matched  to  the  nonlinear 
slowly  varying  oscillations,  resulting  in  the  determination  of  accurate  analytic  formulas  for  the 
boundaries  of  the  basin  of  attraction  and  connection  formulas  across  the  separatrix  for  the  amplitude 
and  phase.  Under  current  investigation  are  generalizations  of  the  slow  crossing  of  a  separatrix  to 
arbitrary  Hamiltonian  systems  and  to  nonchaotic  situations  in  which  small  periodic  forcing  causes  the 
existence  of  an  infinite  sequence  of  resonance  layers  that  coalesce  on  the  separatrix. 


2.  Status  of  Research 


During  the  period  of  the  first  NSF  grant  (1985-87),  research  on  two  problems  was  completed: 

1.  The  leading  tail  of  solitary  waves  for  the  Korteweg-de  Vries  equation  [10] 

2.  Structure  of  two-dimensional  diffusive  shock  waves  [1] 

In  addition,  preliminary  work  began  on  two  problems: 

1.  Crossing  of  a  separatrix  (see  subsequent  discussion) 

2.  Nonlinear  penumbral  caustics  [7] 

We  will  present  a  brief  summary  of  our  results  since  1987  on  the  following  two  related  problems: 

1.  Phase  shift  for  perturbed  strongly  nonlinear  oscillators  and  dispersive  waves 
(summarized  in  greater  detail  in  [15]) 

2.  Analytic  formulas  for  capture: 

the  boundaries  of  the  basin  of  attraction  and  the  slow  transition  across  a  separatrix 
(summarized  in  greater  detail  in  [13]). 

Other  completed  work  concerns  slowly  varying  traveling  waves  for  reaction-diffusion  equations  [16]. 

Major  efforts  continue  on  the  crossing  of  a  separatrix.  An  investigation  ([A]  with  £.  Ho)  has  begun  on 
general  Hamiltonian  systems  under  quasi-autonomous  dissipative  perturbations.  So  far  it  has  been 
shown  that  all  of  the  work  of  Bourland  and  Haberman  on  separatrix  crossing  can  be  extended  to  the 
more  general  case.  Formulas  for  the  phase  shift  have  been  determined.  Accurate  criteria  for  capture 
has  been  obtained,  that  is,  similar  analytic  formulas  exist  for  the  perturbed  boundaries  of  the  basin  of 
attraction.  Numerical  computations  have  motivated  an  asymptotic  investigation  of  higher  order 
corrections  to  the  usual  Melnikov  energy  bounds.  Preliminary  work  has  also  been  done  on  the 
Hamiltonian  system  that  describes  the  usual  weakly  nonlinear  resonance.  An  additional  investigation 
([B]  with  J.  Brothers)  has  begun  on  the  effect  of  small  periodic  perturbations  on  strongly  nonlinear 
oscillatory  dynamical  systems  with  double-well  potentials.  The  object  is  to  extend  the  ideas  of 
Bourland  and  Haberman  for  the  crossing  of  a  separatrix  by  including  a  small  but  fast  periodic  forcing. 
Chaos  can  result  (with  a  corresponding  fractal  boundary  of  the  basin  of  attraction),  but  our  interest 
will  be  on  the  nonchaotic  behavior  preliminary  to  the  larger  amplitude  bifurcation  to  chaos.  The  usual 
slow  variation  procedures  fail  at  each  resonance  band  where  the  forcing  frequency  is  a  fractional 
multiple  of  the  nonlinear  natural  frequency.  There  are  an  infinity  number  of  such  subharmonic 
resonance  bands  that  coalesce  on  the  unperturbed  homoclinic  orbit.  We  have  carefully  calculated  the 
strongly  nonlinear  oscillations  away  from  each  resonance  using  multi-phase  averaging.  In  addition  we 
have  calculated  the  solution  in  each  resonance  layer  and  matched  the  solution  across  each  resonance 
layer  in  order  to  calculate  the  phase  shift.  The  forcing  amplitude  acts  as  a  bifurcation  parameter.  For 
sufficiently  small  amplitudes  (the  ones  we  wish  to  study  at  first)  subharmonic  resonant  periodic 
solutions  do  not  exist  and  the  solution  merely  makes  a  slow  passage  through  resonance  (called  transient 
resonance).  The  more  frequently  studied  case  of  sustained  resonance  where  periodic  solutions  udsts  will 
not  be  particularly  harder  to  analyze,  but  we  still  restrict  our  attention  at  first  to  the  more  elementary 
case.  A  separate  analysis  in  the  neighborhood  of  the  unstable  periodic  solution  will  be  necessary  but 
will  be  quite  similar  to  the  analysis  of  the  homoclinic  orbit  or  separatrix  of  an  unstah  .e  fixed  point  as 
performed  by  Cary,  Escande,  and  Tennyson  and  looked  at  in  more  detail  by  Bourland  and  Haberman. 
We  are  planning  to  do  this  type  of  analysis  so  that  the  near  homoclinic  solution  can  be  matched  to  the 
slowly  varying  solution  away  from  the  homoclinic  orbit.  Again,  great  care  must  be  entertained  since 
the  subharmonic  resonance  bands  coalesce  on  this  homoclinic  orbit.  We  hope  to  obtain  formulas  for 
the  boundaries  of  the  basin  of  attraction  in  the  case  in  which  the  solution  is  not  chaotic.  We  anticipate 
that  a  proposal  to  continue  work  on  this  project  will  be  forthcoming  in  early  Fall  1992. 


3.  Phase  Shift,  for  Perturbed  Strongly  Nonlinear  Oscillators 
and  Dispersive  Waves 

1.  Method  of  multiple  scales.  Modulations  of  traveling  waves  for  nonlinear  partial  differential 
equations  can  be  obtained  by  the  method  of  multiple  scales  in  which  the  solution  is  assumed  to  depend 

T'N 

on  a  fast  phase  ip  =  — -  +  <£(X,T)  and  the  slow  variables  X  =  c  x  and  T  =  e  t,  where  ^(X,T)  is 
the  phase  shift.  To  leading-order,  a  nonlinear  ordinary  differential  equation  in  the  traveling  wave 
coordinate  ip  is  obtained,  from  which  one  derives  the  local  amplitude  dependent  dispersion  relation.  At 
0(e)  in  the  perturbation  expansion,  L(u^)  =  Rp  where  L  is  the  linearized  operator  and  Rj  is  the 
corresponding  nonhomogeneous  term.  Eliminating  the  secular  terms  yields  the  leading-order  wave 
action  equation.  It  is  more  difficult  to  eliminate  secular  terms  from  the  O(e^)  equation,  L(u2)  =  R^, 
but  in  [3]  it  was  shown  that  the  resulting  solvability  condition  determines  the  phase  shift.  It  was  also 
shown  there  that  an  easier  way  to  derive  the  phase  shift  equation  is  based  on  Whitham’s  exact 
equation  for  wave  action: 

FT  1  +  W  ^  =  ~  D’ 

where  I  is  the  exact  wave  action,  Q  its  flux,  and  D  its  dissipation.  If  a  perturbation  expansion  is 
introduced  (I  =  Iq-NIj+...,  Q  =  Qg+«Qj+...,  and  D  =  Dg+eD^-f ...),  then  the  amplitude  parameter  E 
is  determined  by  Ig  +  ^  Qg  =  —  Dg.  Evaluating  the  exact  wave  action  equation  to  0(e)  yields 

w'i  +  lcQi=  -DP  <L2> 

Ip  Qp  and  Dj  depend  on  u^  (the  higher  order  perturbation  of  the  solution).  However,  if  the  O(e) 
perturbation  is  purely  dissipative,  then  in  [3]  it  was  shown  that  0(e)  perturbations  of  wave  action,  its 
flux,  and  its  dissipation  are  only  due  to  perturbations  of  the  wave  number  k  and  frequency  w: 

Ij  =  D(Ig),  Qj  =  D(Qg),  and  Dj  =  D(Dg)+  J g  u^^dip,  (1-3) 

where  D  =  ^  +  kj  ^  ^  is  the  linearized  operator  corresponding  to  a  Taylor 

1  o 

series  in  u  and  k.  The  term  J  g  Ug  d^>  represents  the  dissipation  due  to  the  O(e^)  perturbation. 

0 

Thus,  the  higher-order  wave  action  equation  yields  a  linear  PDE  for  the  phase  shift: 


dn  K  dT  ‘Ow+  dX 


V+M(-^qow+m  Qok 


Qo.)=  -(- 


X 

<‘-4» 


generalizing  earlier  work  on  ordinary  differential  equations  ([2]  and  [6]).  For  non-purely  dissipative 
perturbations,  the  governing  equation  must  be  slightly  modified. 


2.  Lagrangian  formulation.  In  [5]  it  was  also  shown  that  the  phase  shift  can  be  determined  this  way 
for  perturbed  nonlinear  partial  differential  equations  formulated  in  terms  of  a  Lagrangian: 

&h  +  &h-h=-<*-'2*’  <2-» 

where  the  Lagrangian  satisfies  L  =  L  (upUx,u;X,T)  and,  for  example,  Lj  stands  for  the  partial 


derivative  of  L  with  respect  to  the  first  argument  ut.  Whitham  showed  that  the  modulation  equations 
correspond  to  the  exact  conservation  of  wave  action  modified  by  the  inclusion  of  dissipation: 

(2.2) 


0  t  _l  JL  t  —  _  n 

3T  w+9x  k  ~  u’ 


1 


where  L  is  the  averaged  Lagrangian  L  =  J  L  dt/>.  Leading-order  wave  action  determines  the  amplitude 

Q  O  0 

in  the  well-known  way(  —  =  —  Dq  ).  For  purely  dissipative  0(e)  perturbations, 

(1.3)  were  again  shown  to  be  valid  in  [5]  so  that  the  next  order  of  the  wave  action  equation  becomes: 


-M'£)1-n  +  ^6\=-6(D0)-|*u0.,.d^ 


0 


(2.3a) 


where  D  =  —  <p^  — I-  <p-^  .  Thus,  the  phase  shift  <p  satisfies  the  following  linear  partial 

differential  equation  written  in  terms  of  partial  derivatives  of  the  averaged  Lagrangian: 


JL 

dT 


(  ffl  dX  +  dX  '  £T  i'wk+  dx  ^kk)  ~  0T  D0W  dX  °Ok  J  gu0^d^‘  (2‘3b) 


Similarly,  the  phase  shift  was  obtained  in  [5]  for  strongly  nonlinear  multi-phase  oscillatory  waves. 


3.  Korteweg-de  Vries  equation.  In  [4]  these  ideas  were  applied  to  oscillatory  single-phase  solutions  of 
the  Korteweg-de  Vries  (KdV)  equation.  Since  there  are  two  amplitude  parameters,  there  are  two  exact 
action  equations,  each  of  which  can  be  utilized.  The  modulation  equations  for  the  phase  shift  can  be 
obtained,  but  consists  of  a  coupled  system  of  two  linear  partial  differential  equations  in  two  unknowns 
since  only  one  amplitude  parameter  can  be  eliminated  using  the  dispersion  relation. 


4.  Generalization  to  two  fast  scales.  The  usual  method  of  multiple  scales  assumes  the  perturbations 
are  slowly  varying  traveling  waves.  This  puts  a  severe  restriction  on  the  perturbations  of  the  initial 
conditions.  To  overcome  this  restriction,  in  [12]  a  method  of  multiple  scales  was  developed  with  two 
fast  scales,  the  usual  traveling  wave  coordinate  and  time  (ip  =  -  -  +  <£(X,T)  and  t  =  t)  and  the 

usual  two  slow  variables  (T  =  e  t  and  X  =  e  x  ).  It  has  been  shown  that  the  partial  differential 
equation  (2.1)  is  equivalent  to  an  exact  equation  of  wave  action: 

+  =  _D-  (4J) 
where  the  Lagrangian  is  now  averaged  over  time  as  well  as  the  phase: 

t  t  1 

I  =  lim  |  f  L  dt  =  lim  if  f  L  dtp  dt;  (4.2) 

t— >oo  1  J  t— >00  1  J  J 

0  0  0 


it  is  the  time  average  of  the  usual  averaged  Lagrangian.  When  the  wave  action  equation  is  evaluated 
to  leading  order,  the  usual  equation  for  the  wave  amplitude  results.  To  evaluate  the  wave  action  to  the 


next  order  requires  knowledge  of  the  perturbation  of  the  solution  u^.  Here,  u^  =  u  j(0)  +  Uj(0,t).  In 
general  the  time  dependence  is  difficult  to  obtain  explicitly  since  it  corresponds  to  the  equations  that 
arise  in  analyzing  the  stability  of  the  traveling  wave.  However,  for  the  nonlinear  Klein-Gordon 
equation,  the  perturbed  wave  action  Ij  satisfies 

t  1 

II  =  D  (IQ)  +  Urn^  {  |  |  (ufl  flj  -  2u>0  ufl  uj  )  dtf  dt.  (4.3) 

0  0  V 

In  [12]  we  showed  the  limit  in  (4.3)  vanishes  due  to  the  boundedness  of  fij  (corresponding  to  the 
modest  assumption  that  the  traveling  wave  is  stable).  Thus, 

h=  6  (IQ)  ;  (4.4) 

the  time-dependent  part  of  u^  does  not  contribute  to  perturbations  of  the  wave  action,  its  flux,  and  its 
dissipation.  Consequently,  (1.4)  or(2.3)  is  not  altered  when  initial  conditions  are  properly  analyzed. 

5.  Method  of  averaging.  For  ordinary  differential  equations,  methods  of  averaging  (based  on  energy- 
angle  variables  and  near  .  ientity  transformations)  have  proved  useful  to  derive  the  long  time  behavior 
of  perturbed  periodic  systems.  Recently,  a  similar  formalism  has  been  introduced  in  [14]  for  oscillatory 
dispersive  wave  solutions  of  nonlinear  partial  differential  equations.  Amplitude  parameter  E  and  angle 
ip  variables  are  introduced,  motivated  by  the  corresponding  expressions  derived  by  the  method  of 
multiple  scales.  Two  “energy”  equations  have  been  derived  in  terms  of  the  appropriately  defined  wave 
number  k  and  frequency  w,  which  are  referred  to  as  standard  form  for  the  amplitude  equations: 

A  [S],  +  it-],  =  ‘l)  +  '2l\  <5-D 

where  A  and  B  are  related  two-by-two  matrices.  The  expressions  for  these  matrices  and  the 
perturbation  terms  can  be  found  in  [14],  Equations  (5.1)  form  a  hyperbolic  system  (which  can  be 
shown  to  have  characteristic  velocities  ±c  as  does  the  pde).  It  has  been  shown  that  the  traveling  wave 
assumption  corresponds  to  u^=  —  «u^  and  ux=ku^,  from  which  the  standard  form  for  the  angle 
equations  has  been  derived: 

rPt=  ~w  +  f0(kt,wt,kx,wx)  +  dl+<2  f2  (5.2a) 

^x  =  k  +  go(kt>wt’kx>wx)  +  €8f  +  <2  «2-  (5<2b) 

Detail  expressions  for  fg  and  gg  also  appear  in  [14].  The  four  equations,  (5.1)  and  (5.2),  which  involve 

the  three  unknowns  (k,w,i/>)  are  linear  in  spatial  and  temporal  derivatives  of  k  and  w.  They  are  well- 
posed  since  the  two  consistency  equations  have  been  shown  to  be  equivalent.  Although  these  equations 
appear  complex,  after  the  introduction  of  near  identity  transformations  and  averaging,  modulation 
equations  for  the  amplitude,  wave  number,  frequency,  and  phase  shift  have  been  derived  [14]  with  some 
effort  which  are  equivalent  to  the  ones  obtained  by  the  method  of  multiple  scales. 


4.  Analytic  Formulas  for  Capture:  the  Boundaries  of  the  Basin  of 
Attraction  and  the  Slow  Transition  Across  a  Separatrix 

1.  Before  capture:  strongly  nonlinear  slowly  varying  oscillations.  To  determine  the  slow  variation  of  a 
nonlinear  oscillator  due  to  small  perturbations,  various  techniques  have  been  developed.  It  is  well- 
known  that  these  asymptotic  techniques  fail  if  the  trajectory  approaches  an  unperturbed  homoclinic 
orbit,  the  infinite  period  limit  of  a  periodic  solution.  For  Hamiltonian  systems  (without  dissipation), 
Timofeev,  Neishtadt,  and  Tennyson,  Cary,  and  Escande  represented  the  solution  in  the  neighborhood  of 
the  unperturbed  homoclinic  orbit  as  a  large  sequence  of  nearly  homoclinic  orbits.  They  showed  that 
the  action  changed  by  a  small,  but  important  amount  after  the  solution  crossed  the  separatrix.  We 
have  analyzed  the  slow  crossing  of  a  separatrix  when  small  dissipation  is  present: 

§  +  Vy(y)  =  -  A*(y,  ,  (1.1) 

where  V  is  a  double-well  potential  with  a  saddle  point  at  y=0.  We  assume  the  perturbation  is  purely 
dissipative,  h*(y,  —  ^)  =  —  h*(y,  ^).  With  the  inclusion  of  small  dissipation,  it  is  well  known 
that  the  two  branches  of  the  stable  manifold  are  tightly  wound  (separating  solutions  captured  in  the 
left  well  from  those  captured  in  the  right).  In  [8]  we  analyzed  the  slow  passage  through  a  separatrix 
due  to  dissipation  using  the  method  of  multiple  scales,  determining  the  boundaries  of  the  basins  of 
attraction  and  connection  formula.  Generalizations  to  slowly  varying  potentials  V(y,tt)  and  dissipative 
perturbations  h*(y,dy/dt,ct)  are  discussed  in  [9].  In  [11]  we  used  the  equivalent  method  of  averaging 
with  the  energy  E  and  angle  ip  variables  (easily  related  to  the  usual  action-angle  variables): 


E=  i(fe2  +  V(y) 


dy 

v/2[E-V(y)]1/2 


(1.2a, b) 


where  fl(E)  is  the  usual  frequency.  The  following  averaged  equations  were  derived  in  [11]  using  a  near 
identity  transformation,  if  the  perturbation  is  purely  dissipative: 


%  =  -e  11(e)  D(e)  +  0(c3)  and  ^  =  fi(e)  +  0(c2), 


(1.3a, b) 


where  e  is  the  average  of  E.  We  introduce  a  perturbation  expansion,  e=eg+€e^+  O(e^).  The  phase  has 
been  derived  in  [2]  using  the  method  of  multiple  scales  and  in  [11]  using  the  method  of  averaging: 


<P  =  \  l  dT  +  A  n0(T)  +  B  +  0(c)  , 
where  Oq  =  n(eg).  From  a  careful  analysis  of  the  initial  conditions  ([2]  and  [11]), 


(1.4a) 


D0(e0(0)) 


hyoti’d^'  +  initial  perturbation  of  the  physical  energy  E)  .  (1.4b) 


The  slow  variation  of  the  leading-order  energy  (equivalent  to  the  usual  action  equation), 

^>  =  -n0D0(e0),  (1.5) 

follows  for  (1.3a).  As  the  homoclinic  orbit  is  approached  (eg  — *  0),  1/Q0  —  Cj  In  |  e0  |  +  c2  ,  where  c2 
is  shown  in  [8]  to  depend  on  the  entire  potential.  For  captured  oscillations  in  the  right  or  left  well,  cj1 
=—  j  2A,  where  A  is  given  by  (3.1),  and  c2  is  different  in  each  case.  The  time  Tc  (when  e0  =  0) 
may  be  determined  from  (1.5).  At  this  time,  the  corresponding  phase  follows  from  (1.4a). 

2.  Capture  and  the  boundary  of  the  basin  of  attraction.  From  (1.1),  dE/dt  =  —  ch*  dy/dt.  As  the 
solution  approaches  the  critical  energy  (e0  — *  0),  the  solution  is  composed  of  nearly  solitary  pulses.  The 
Melnikov  result  is  that  the  energy  decreases  by  eD^  >  0  (eD^  >  0)  over  a  right  (left)  homoclinic  while 
over  a  complete  (figure-eight  shaped)  pulse,  it  decreases  by  f  •  The  last  saddle  approach 

(see  §  3)  has  positive  energy,  W0  >  0;  the  energy  will  become  negative  at  the  next  saddle  approach. 
Since  the  energy  will  diminish  by  cDR  on  the  next  right  homoclinic  orbit  (the  well-known  separation  of 
the  two  stable  manifolds  near  the  separatrix),  the  criteria  for  right  capture  is  that  0  <  W0  <  cDft. 
Corresponding  initial  conditions  can  be  determined  from  slow  variation  theory,  since  when  the  phase  is 
an  integer  we  ([8],  and  [11])  have  shown  that  eg  of  the  method  of  averaging  and  the  actual  energy  E 
are  sufficiently  close  near  the  separatrix.  (At  other  phases,  the  difference  is  O(e),  which  is  too  large.) 
Near  the  separatrix  for  the  double-pulsed  solution,  the  dissipation  Dg  — ►  Dp+D^,  so  that,  from  (1.5), 
Aeg/A0  =  —  (Dp+D^).  We  use  the  time  of  the  last  saddle  approach  (where  the  phase  is  an  integer) 
and  the  time  Tc  (where  eg  =0  and  the  phase  equals  <f>c).  When  the  phase  is  an  integer,  eg  is 
sufficiently  close  to  the  actual  energy,  which  is  Wg  +  if  using  the  Melnikov  idea.  Thus, 

W0  +  ^dl=  <(dr+  D^C*,  (2>1) 

where  =  <j>c  -  [<£J  and  where  [<t> J  is  the  integral  part  of  <j>c.  By  an  appropriate  integration  of 

(1.5),  we  determine  the  initial  energy  e(0)  whose  energy  at  the  last  saddle  approach  is  Wg: 

1  f  (  \  dE  -  U]+  IPL+  (Wq/c)  _  1  *h  dl()  (22) 

7  J  D0(E)  “  ^  +  Dr  +  Dl  Do(e(0))  J  hy°*  ^  '  (  ’  } 


hy0w,di(-  , 


where  (here  only)  the  perturbation  of  the  initial  energy  (see(1.4b)  is  assumed  to  vanish.  The 
boundaries  of  the  basin  of  attraction  (for  the  right  well)  correspond  to  Wg  =  0  and  Wg  =  c  DR. 

3.  Separatrix  Region:  a  Sequence  of  Solitary  Pulses.  A  transition  region  near  the  unperturbed 

homoclinic  orbit  consists  of  a  large  sequence  of  solitary  pulses  (separated  by  a  close  approach  to  a 

saddle  point),  as  discussed  for  Hamiltonian  systems  (without  dissipation)  by  Cary,  Escande,  and 

Tennyson.  The  solution  in  the  vicinity  of  a  saddle  point  with  energy  Wn  is 
IW  if  , 

y  =  ±  - — r* —  sinh  or  cosh  A>/2(t-tn),  where  A2=~V  (0)  >  0, 


(3.1) 


and  where  t  =  tQ  is  the  time  of  closest  approach  to  the  saddle.  The  energy  dissipation  depends  on 
whether  the  solution  is  nearly  a  right  or  left  homoclinic  orbit,  Wn_j  — Wn  =eDR  or  eDL  .  To  determine 
the  time  of  the  saddle  approaches,  we  consider  the  approximate  homoclinic  orbit  connecting  two  saddle 
approaches  occurring  at  t=tQ  _  ^  and  t=tn.  The  exponential  growth  specified  by  (3.1)  must  correspond 
to  the  exponential  decay  (calculated  in  [8])  of  the  tails  of  the  corresponding  solitary  pulses.  In  this 
manner,  the  “period”  for  a  sequence  of  solitary  pulses  is  the  average  of  the  periods  for  the  surrounding 
saddles: 

J2  1 /2 

,n-t»-l  =  -Trl"lWnWn-ll  +  *1  •  <3-2)  - 

where  for  left  homoclinic  orbits  c2= c2  in  (3.2)  .while  for  right  ones  c2=  c^\  It  is  not  difficult  to  solve 
(3.2)  in  terms  of  the  unknown  t0  using  the  basic  property  of  the  gamma  function,  T(x  +  1)  =  xr(x). 

4.  Matching  and  the  Connection  Formulas  Across  a  Separatrix.  After  the  solution  is  captured  into  the 
right  well,  we  can  insist  that  eQ  be  the  unique  solution  of  (1.5)  which  starts  at  the  critical  energy, 
Cq(Tr)  =  0,  at  some  as  yet  unknown  time  T=TR=  ftR.  The  phase  (see  (1.4a)  contains  the  second 
constant  0(TR).  Using  the  same  reasoning  used  to  obtain  (2.1),  we  derive  <fi(TR)=W0/e  DR. 
Eliminating  W0  from  (2.1),  yields  DR<^(TR)  =  ^  (DR+  DL)  —  ^  DL,  determining  the  captured  phase 
ftom  initial  conditions.  The  asymptotic  expansion  for  nonlinear  oscillators  fails  when  e^  =0(  ),  i.e.  as 
the  unperturbed  separatrix  is  approached,  suggesting  the  method  of  matched  asymptotic  expansions 
(even  though  the  energy  is  predicted  correctly  at  extrema).  Similarly,  the  sequence  of  solitary  pulses 
will  fail  when  n=0(  l/«),  since  there  the  solution  can  no  longer  be  ipproximated  by  homoclinic  orbits. 
We  [8]  have  shown  that  the  energy  already  matches  to  leading  order  before  capture  due  to  (2.1),  which 
determines  Wq  .  The  careful  matching  [8]  of  the  phase  angles  before  capture  determines  the  time  t0  of 
the  last  saddle  approach.  In  the  same  manner,  the  phase  angle  near  homoclinic  orbits  was  matched  to 
the  phase  angle  of  the  nonlinear  oscillator  after  capture,  determining  the  starting  time  Tj^  =  e  tj^  : 


where  Cj  =  c^+ci  =2c^=— \/2/A  and  Dq=Dj^+D^,  completing  the  connection  formulas  for  the  slow 
passage  through  the  separatrix. 
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Hamiltonian  systems,  (b)  analytic  formulas  for  the  boundary  of  the  basin  of  attraction  for  Hamiltonian 
systems,  (c)  capture  in  1-1  resonance,  (d)  logarithmic  correction  to  Melnikov  energy  dissipation  for 
autonomous  systems,  and  (e)  logarithmic  correction  to  Melnikov  energy  dissipation  using  time- 
dependent  methods. 


It  is  anticipated  that  R.  Haberman  will  write  papers  with  J.  Brothers  on  (a)  multiphase  averaging,  (b) 
phase  shift  jumps  across  sequences  of  resonant  layers,  and  (c)  crossing  a  separatrix  with  periodic 
forcing. 
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